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MULTIPLE  SCATTERING  EFFECTS  IN  RADAR  OBSERVATIONS  OF  WAKES 


I.  Introduction 

%  I 

The  large  amplitude  waves  observed  In  ship  wakes  and  the 
large  radar  returns  Including  sharp  angular  dependences  suggest  that 
first  order  Bragg  scattering  theory  Is  inadequate  to  describe  the 
experiments. 

Here  we  consider  a  simple  theory  to  qualitatively  investigate 
higher  order  effects.  There  are  two  important  limitations.  . 

(a)  We  restrict  our  consideration  to  a  simple  scalar  problem 
with  simple  boundary  conditions.  A  full  scale  theoretical  treatment 
of  the  true  electromagnetic  problem  can  readily  be  carried  out. 
However,  it  is  much  more  complicated,  probably  obscures  the  central 
points,  and  in  view  of  (b)  probably  unwarranted  at  this  point. 

(b)  The  detailed  nature  of  the  scattering  wave  field  is  very 
poorly  known.  To  get  some  insight  we  have  applied  the  formula  to  an 


Idealized  mathematical  description  of  a  Kelvin  wake  and  to  a  model 
"derived1*  from  the  Dabob  Bay  experiments. 


The  essential  conclusions  are  that  there  can  be  look  angles 
such  that  second  order  scattering  can  be  larger  than  first  order  and 
comparable  in  magnitude  with  the  first  order  scattering  When  that  is 
significant. 

tfe  emphasize  the  weakness  of  the  conclusions  because  of  the 
lack  of  detailed  knowledge  of  the  basic  wave  fields.  Recommendations 
to  obtain  the  needed  information  are  made. 
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A  Staple  Theory 


He  consider  the  electromagnetic  field  to  be  described  by  a 
simple  scalar  field  satisfying  the  wave  equation. 


(V2  +  k2)  +  -  0  . 


(1) 


An  incident  field 


♦i  “  e 


i(k  •  r  -  k  Z) 


2  2  2 
<k^  +  k  )  -  k 
~h  v 


(2) 


Is  to  be  scattered  by  a  patch  of  a  wavy  surface  centered  on  the  plan 
Z  “  0,  l.e.,  the  surface  S  Is 


Z  -  h(x,y)  (3) 


where  h  Is  to  be  small. 


As  a  simple  boundary  condition  we  choose  $  *  0  on  S  . 


We  decompose  the  total  field  $  Into 


♦  ■  ♦jj  +  A$ 


(4) 


Here  ♦  is  the  field  which  would  result  if  h  ■  0. 
o 


r  -i  k  Z  i  k  Z 


i  lc.  •  r  _ 

,  .  ~h  ~  ,  v  v  . 

i.e.  ♦  -  e  {e  -e  } 


(5] 


Further  we  introduce  a  Green's  function  G  (r'»  r)  satisfying 


,.2 


^  G  (r',  r)  -  6(r  -  r') 


(6] 


and  G  (x',  y',  o;  r  )  «  0 


(7] 


Applying  Green's  theorem  we  obtain 


(r)  -  s/  {  A*(r'>  JL  G(r',r)  -  G(r',£)  |£±  (r') 


3n' 


dS' 


(s: 


For  a  perturbation  theory  we  assume  h  -  e  and  consider  a 
formal  expansion  in  c  • 


Then  A$(r)  -  e^Cr)  +  c  ^(r)  +  .  . 


(9 


Restricting  ourselves  up  to  second  order  terms  this  can  be 


mitten 


A*(r)  - 


But 


and 


/  /  dx'dy'  A*<0  G  (r'.r)  -  G(r',r)  ||* 
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We  readily  construct 


i  [K  (x-x')  +  K  (y-y') 

_  ,  .  %  -i  ,  d  K  e  7 

G  (r  ,r)  -  - -  f  - — -  l  > 


(2*)‘ 


2u 


ip(z-z')  -iii(z-z')  iii(z+z")  -iu(z+z') 
{  }  »  e  +e  -e  -e 


Then  to  first  order  we  have 


^(r)  -  //  dx'dy'  ^(x'.y'.O) 

e^Cx',  y',  h)  +  ^(x',  y',  h)  -  0 

ik  •  r' 

^(x',  y',  0)  -  2ikyh  e 

■—  G(x',  y',  o;  r) 

,  i  [k  (x-x')  +  k  (y-y')] 
-=  fl  e  x  7 


G(x',y',0;r) 
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(12 
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. .  ,,,,,,  1  r  .2  ,  i(V  +  kyy)r  iuZ  -iuZ, 

♦  .(r)  -  2ik  SI  dx  dy  - =-  Jdke  J  [e  ^  -  e  ] 

V  (2*/ 

i[ (k1  -  k  )  x'  +  (k1  -  k  )  y'] 
x  h(x',  y')  e  x  y  y  y 


Now  if  we  write 


1  ,  ,2 


h(x',  y') - i-r  /  d  q  h  (a'>  e 


_i  a'*  f 


1  2  .  •  IJ  4»»<va>z  -^(Va>z  , 

♦  (r)  -  2ik  — dVh(a')  e  -1  (e  -e  ) 

V  (2w)Z 

(15) 


From  Equation  (8)  we  then  read  off  that 


♦,(r)  -  IS  dx' dy'  Mx',  y',  0)  - 5- G  (x',  y',  0,  r) 

z  3Z 


+  4>,  (x',  y'»  0)  -%pr  G(x',  y',  0;  r) 


-  h  G(x' ,  y',  0;  r)  7—  ^  (x',  y',  0) 


Note:  With  our  boundary  condition 


Also 


0 


3  G  (x',  y't  0;  r)  -  0  . 

az'2 

e+j^Cx'.  y',  h)  +  e2  ^(x1,  yl,  0)  +  ^  - 
2 

But  to  order  e  we  have  - |  ■  0 

3Z"4  Z'  -  0 

♦2(*'»  y'*  °>  — h  *1  <*'•  y'*  0) 

and  so 

*2  -  -2  II  dx'dy'  h  *1(x',  7' ,  0)  —  G(x',  y',  0;  r)  (16) 

which  yields 

~®*kv  2  2~  ~  D 
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<2*)4  x 

ii)5i  -  a  -  el  •  £  ,  1p<k1  -  a  -  k>z  _1p(ki  -  a  -  k>z 

x  e  I  e  -e 


In  particular  If  we  ask  for  the  amplitude  of  the  wave  In  the 
back  scattered  direction,  we  obtain 


(a)  From  equation  (15) 


III.  Application  to  "Mathematical"  Kelvin  Wakes 

Let  us  now  apply  these  results  to  the  Kelvin  wake  as  computed 
by  simple  stationary  phase  methods.  There  are  many  things  wrong  with 
this  simplest  picture  of  the  wake,  but  we  believe  that  it  gives  a  good 
enough  picture  of  the  wave  heights  for  an  analysis  of  the  importance 
of  multiple  scattering. 


The  stationary  phase  picture  asserts  that  at  a  point  (x,y)  in 
the  wake,  making  an  angle  6  with  respect  to  the  wake  axis,  the 
dominant  wave  vectors  in  the  surface  height,  h,  are 


where  U  is  the  ship  velocity  which  is  in  the  x  direction.  This  means 
that  the  ocean  wave  height  function  is 


h(x,y)  =*  A(x,y)  cos  4>  (x,y) 


<f(x,y)  =  k  (x,y)  x  +  k  (x,y)  •  y  . 

•*  y 
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where  A(x,y)  Is  a  slowly  varying  function,  poorly  calculable  from 
stationary  phase  arguments,  and  4  Is  rapidly  varying  and,  of  course, 
determined  by  the  stationary  phase  arguments.  There  are  actually  two 
stationary  phase  solutions  for  &  :  the  one  given,  corresponding  to 
the  diverging  wave  train,  and  another  one  corresponding  to  the 
transverse  wave  train.  The  two  trains  are  displayed  In  Figure  1.  The 
short  wavelengths  of  direct  Interest  In  radar  backscatter  are  to  be 
found  In  the  diverging  wave  system,  so  we  will  not  concern  ourselves 
here  with  the  transverse  waves.  We  will  later  resort  to  a  combination 
of  theoretical  argument  and  direct  observation  to  determine  reasonable 
values  for  A  In  various  parts  of  the  wake. 

To  model  a  SAR  observation  of  the  ocean  surface,  we  assume 
that  the  radar  processing  In  effect  forms  a  beam  which  Illuminates  a 
patch  on  the  surface  of  linear  dimension  b  centered  on  the  point 
(x  ,y0).  The  radar  return  Is  therefore  computed  from  the  formula  of 
the  previous  section  by  extending  the  spatial  Integrals  only  over  the 
window  x  -  y  <  x  <  x  +y,y  y  <  y  <  y  +  y  •  The  scattered 
intensity  will,  of  course,  depend  on  the  patch  size,  b. 

The  scattering  formula  are  expressed  in  terms  of  the  Fourier 
transform,  h(q)  ,  of  the  ocean  surface.  By  the  preceding  remark,  we 
want  to  compute  the  transform  of  the  wave  height  times  the  radar 


Figure  1 .  Kelvin  Waves 
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window  function.  The  result  for  a  patch  centered  on  a  point  at  wake 
angle  6  is 


h(q>  -  A(9)  dedn  cos  (-  (-8=-  -  q  )e  +  <-fy  -  o)n] 

D/Z  2D  0  4 U  0  T 

In  writing  this  expression  we  have  simplified  the  variation  of  (k  ,k  ) 

x  y 

across  the  radar  patch.  The  omitted  effects  would  contribute  terms 
quadratic  in  e  and  n  to  the  phase  of  the  cosine.  We  have  verified 
that  for  parameter  values  appropriate  to  the  SEASAT  or  Dabob  Bay 
problems,  this  approximation  is  reasonable.  The  result  is 

h(q?8)  -  A(8)b2  sine  [(q  +  -&y)  y]  sine  [(q - §-y)  y 

*  2V  d  Z  y  4U  6 

For  our  purposes,  it  suffices  to  know  that  h(q,8)  has  a  maximum 
2 

amplitude  of  b  A(9)  and  that  it  is  peaked  at  the  appropriate 
stationary  phase  wave  number  with  a  width  in  wave  number  space  of 
2ir/b  . 


We  are  now  able  to  compute  the  radar  backscatter  amplitude  in 
first  and  second  Born  approximation.  According  to  the  previous 
section, 


-  C  h  (2kh,9) 
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-1C  / 


4  /  k  -  q 


h(yq,  0)  h  (kh  -  q,  0) 


4> 


2 


(2n)2 


where  k  Is  the  magnitude  of  the  radar  wave  vector,  k^  is  the 

horizontal  projection  of  that  wave  vector  and  C  is  a  common  factor  of 
2 

dimension  1/L  which  we  would  need  to  know  in  order  to  get  the 
absolute  scattered  intensity. 


Since  h  is  a  sharply  peeked  function  of  its  argument,  2 
will  be  large  only  for  a  narrowly-defined  band  of  0  .  If  we  define 

*KELVIN(0)  ”  ( - 2~  *  ~~2~2) 

RELVIW  21T0  4UZQ4 

then  the  condition  for  first  and  second  order  scattering  to  be  large 
is 


first  order:  2^  - 

second  order:  ^  -  *^^(0) 

The  angular  width  of  the  first-order  Born  pattern  is  determined  by  the 
width  in  momentum  space  of  h  and  the  rate  of  variation  with  0  of 
^KELVIN*  A  bit  algebra  shows  that 


2* 


50  _ 

9  ’  b  ‘  ^KELVIN 

If  b  ~  6  m  and  ~  «3  m  (typical  values  for  SEASAT),  this 

gives  66/8  ~  1/20  .  This  is  very  narrow  indeed,  but  not  out  of  line 
with  observations. 

Let  us  finally  consider  the  relative  sizes  of  first  and 
second  order  scattered  amplitudes.  Due  account  being  taken  of  the 
peek  amplitude  of  h  and  its  momentum  space  width,  we  find  that  the 
maximum  backscatter  amplitudes  in  first  and  second  Born  approximations 
are 


♦  “  -  c  •  b2A(V 

•j™1  -  (2kA<e2>)  •  C  •  b2A<82) 

8 j  and  6^  are  the  (different  1)  wake  angles  at  which  the  two  types  of 
scattering  reach  their  maximum,  k  is  the  radar  wave  vector  and  b  is 
the  radar  patch  size.  The  integral  over  q  in  the  expression  for 
is,  because  of  the  narrowness  in  wave  vector  space  of  the  function 
h  concentrated  at  q  -  0  ,  leading  to  a  considerable  simplification 
in  the  formulas. 
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We  expect  that  A(8)  is  a  reasonably  slowly  varying  function 


of  8  ,  except  in  the  near  neighborhood  of  8  ■  0  .  Therefore,  to 
compare  first  and  second  order  scattering  we  will  assume  that 
ACS^)  -  A(8£)  ,  obtaining 


4 

— *  4kA  ■  8itA/X 

radar 


From  slender  ship  theory  we  extract  the  result  that  the  amplitudes  of 
the  diverging  and  transverse  wave  trains  are  related  by 


A  ^ 

diverging 


,3/2 


transverse 


for  8  not  too  small.  One  readily  finds  transverse  wave  amplitudes  of 
.5  m  a  kilometer  behind  a  large  ship.  This  would  correspond  to  a 
diverging  wave  amplitude  of  1.5  cm  at  8  •  6°.  We  will,  therefore, 
take  1  cm  as  a  representative  value  for  A(  8  ).  This  is  confirmed  by 
the  in  situ  wave  height  measurements  of  the  Dabob  Bay  experiment.  For 
a  radar  wavelength  of  23.5  cm  (appropriate  to  all  the  measurements 
under  discussion)  we  therefore  have  ~  *  • 


Our  conclusion  is  that  surface  ship  Kelvin  wakes  are  strong 
enough  that  multiple  scattering  from  the  wake  itself  cannot  be 


neglected. 


At  a  minimum,  this  means  that  multiple  returns,  due  to 
various  orders  of  scattering  ought  to  be  seen.  This  might  be  the 
explanation  of  the  multiple  V  structures  seen  in  some  SEASAT  pictures. 


IV 


Dabob  Bay 


In  Che  Dabob  B ay  experiment,  the  radar  return  from  a  wake  was 
measured  and  the  actual  wave  height  was  measured  in  a  one-dimensional 
transverse  cut  across  the  wake.  The  maximum  radar  return  came  from 
points  at  wake  angle  8  ~  3.6°  .  The  direct  measurement  of  the  wave 
height  in  this  region  Indicates  that  there  is  little  energy  at  wave 
numbers  appropriate  to  first  order  scattering  and  a  substantial  peak 
at  wave  numbers  appropriate  to  second  order  scattering.  This  peak 
does  not  occur  at  the  wave  number  predicted  by  simple  Kelvin  theory. 
This  is  a  hydrodynamic  puzzle  about  which  we  can  say  nothing.  The 
question  of  Interest  to  us  is  whether,  given  the  measured  surface 
height,  it  is  reasonable  that  second  order  scattering  should  overwhelm 
first  order  scattering. 

The  r.m.s.  displacement  recorded  at  the  interesting  wake 
angle  appears  to  be  about  1  cm,  and  corresponds  to  surface  wavelengths 
of  order  40  cm.  As  was  noted  in  the  previous  section,  the  quantity 
2kA,  which  determines  the  relative  amplitude  of  second  and  first  Born 
scattering  for  equal  rms  surface  displacement,  is  of  order  1  for 
A  1  ca  and  *ra<jar  “  24  cm.  An  examination  of  the  power  spectral 
density  of  the  surface  slope  shows  that  the  p.s.d.  is  down  by  a  factor 
of  10~2  from  the  peak  value  at  the  wave  number  appropriate  to  first 
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order  Born  scattering.  Consequently,  the  intensity  of  second  order 

o 

Born  scattering  should  be  about  10  times  that  of  first  order  Born 
scattering. 

So,  if  we  modify  our  scattering  theory  to  include  higher 
order  effects,  there  is  no  inconsistency  between  the  location  of  the 
intense  radar  returns  within  the  wake  and  the  directly  observed  wake 
structure.  Obviously,  one  should  combine  the  observed  wake  structure 
with  a  radar  scattering  model  to  obtain  a  predicted  radar  wake 
structure  and  attempt  a  detailed  comparison  with  the  observed  radar 
wake  images.  An  important  missing  element  is  a  two-dimensional  map  of 
the  wake  structure  since  that  is  what  really  enters  the  radar  cross- 
section.  We  have  used  a  one-dimensional  wake  crossing  measurement  for 
rough  comparison  purposes,  but  must  be  aware  that  surprises  may  occur 
when  we  confront  the  two-dimensional  data  derived  from  photographic 
Images  of  the  wake. 
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V. 


Conclusions 


These  are  based  on  a  simplified  scattering  model  and  some 
questionable  assumptions  as  to  the  nature  of  the  scattering  field. 

The  scattering  model  can  be  readily  improved.  Thus  with  an  increase 
in  the  complexity  of  the  formula  intrinsically  electromagnetic  effects 
(such  as  polarization,  Fresnel  diffractions,  and  dielectric  properties 
can  be  Included.  This  hardly  seems  warranted  till  more  is  known  about 
the  hydrodynamic  structure  of  the  wake. 

Two  hydrodynamic  descriptions  have  been  used. 

(a)  The  idealized  mathematical  Kelvin  wake. 

(b)  Some  estimates  have  been  made  using  the  Dabob  Bay  data. 

Our  conclusion  is  that  first  order  Bragg  theory  is  not 
generally  adequate.  At  some  viewing  angles  there  will  be  strong  first 
order  Bragg  scattering  with  a  sharp  angular  dependence.  At  other 
angles  first  order  Bragg  scattering  will  be  small.  However,  there  can 
be  second  order  Bragg  scattering  of  an  intensity  comparable  to  that 
where  first  order  is  large.  The  angular  width  of  the  second  order 
observations  should  be  greater  than  that  of  first  order,  but  not 
necessarily  dramatically  so. 


The  main  limitations  on  our  present  work  is  that  only  special 
directions  of  look  are  adequately  described.  Dabob  Bay  experiments 
are  closest  to  what  are  considered. 

To  compute  what  should  be  observed  at  other  aspect  angles 
empirical  two  dimensional  wave  height  spectra  are  needed.  Given  these 
the  conclusions  dram  about  angular  sharpness  might  be  modified. 

These  also  might  help  us  to  understand  the  experiment  observations  at 
different  aspect  angles. 

Some  of  the  needed  information  could  perhaps  be  obtained  by 
reworking  the  Dabob  Bay  data.  Also  some  of  this  might  be  obtained 


from  the  NOSC  data 
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